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Abstract 

In previous articles we have developed a theory of down conversion in 
nonlinear crystals, based on the Wigner representation of the radiation field. 
Taking advantage of the fact that the Wigner function is always positive in 
parametric down conversion experiments, we construct a local hidden vari- 
^ | ables model where the amplitudes of the field modes are taken as random 

variables whose probability distribution is the Wigner function. In order to 
achieve our goal we give a model of detection which is fully local but departs 
from quantum theory. In our model the zeropoint (vacuum) level of radiation 
lies below a threshold of the detectors and only signals above the threshold 
are detectable. The predictions of the model agree with those of quantum 
mechanics if the signal intensities surpase some level and the efficiency is low. 
This is consistent with the known fact that quantum mechanics is compatible 
with local realism in that case (a fact called the "efficiency loophole"). Our 
model gives a number of constraints which do not follow from the quantum 
theory of detection and are experimentally testable. 
PACS number(s): 42.50.Ar, 03.65. Sq, 42.50.Lc 



1 



1 Introduction 



The aim of the present article is to propose a local hidden variables (LHV) 
model for the experiments involving parametric down conversion (PDC). The 
model is based on the quantum theory of PDC when formulated in the Wigner 
representation. The reader may wonder what is the use of such a LHV model 
and therefore this is the first question which we should answer. Before doing 
that, we shall briefly comment on the present status of the empirical tests of 
LHV theories. 

As is well known, in 1964 Bell proved that there are predictions of quan- 
tum mechanics which cannot be reproduced by any LHV theory J]]. There- 
fore it seems possible, in principle, to discriminate empirically between quan- 
tum mechanics and "local realism" (i.e. the whole family of LHV theories) 
by means of some specific experiments. Such experiments usually attempt to 
test whether a Bell's inequality is violated and will be called "Bell's tests" in 
the following. Actually most of the Bell's tests made during the last 20 years 
have used photon pairs produced in the process of PDC @, 0, f§. PDC has 
also become popular for the study of other nonclassical aspects of light || || 
and, more recently, for the development of the quantum theory of informa- 
tion like the implementation of cryptographic schemes J7|, teleportation [§], 
etc.. In general the performed Bell's tests have confirmed quantum mechan- 
ics but, in spite of the great effort made, they have been unable to provide 
an uncontroversial disproof of local realism. The reason is that, up to now, 
nobody has been able to test empirically a genuine Bell's inequality, i.e. an 
inequality derived from the assumptions of realism and locality alone. All 
inequalities actually tested in the performed experiments involve additional 
assumptions, allegedly plausible, like fair sampling, no- enhancement, etc. [§], 
|T0|j . The consequence is that the important question whether LHV theories 
are possible is still open. This simple fact is rarely acknowledged in the cur- 
rent literature. On the contrary, it has been repeatedly claimed, in respected 
books and journals over the last thirty years, that local realism has been 
refuted experimentally. 

The essential difficulty for the performance of genuine Bell's tests with op- 
tical photons has been the low efficiency of photon counters, which allows for 
a loophole in the disproof of LHV theories known as the "efficiency loophole" . 
During the eighties and early nineties there was the hope that the efficiency 
loophole could be easily blocked when detectors of high efficiency (and low 
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dark rate) were developed. In the last few years such detectors have become 
available but nevertheless all attempts at blocking the efficiency loophole 
have failed |TT[ and people are turning to experiments not involving optical 



photons |T!| . This problem leads us naturally to ask why it is so difficult 
to block the efficiency loophole in PDC experiments. One of the purposes 
of the present article is to provide a partial resolution of this mystery. In 
fact, we shall prove that it is possible to find a LHV model for all PDC 
experiments by combining standard quantum theory for the production and 
propagation of light with a non-quantum model of detection. Our detection 
model, although not quantal, agrees with quantum theory for the performed 
experiments (which, as said above, suffer from the detection loophole) but 
departs from quantum mechanics in some (as yet unperformed) experiments, 
which might therefore be able to discriminate between quantum mechanics 
and local realism. In general, a LHV model for an experiment, or a class 
of experiments, has two goals. Firstly it proves that LHV theories for those 
experiments do actually exist, a proof which cannot be derived just from the 
fulfillment of Bell's inequality because these are necessary, but not sufficient, 
conditions for local realism. Secondly LHV models play the role of coun- 
terexamples. That is, their mere existence proves that the said experiment is 
unreliable for the disproof of the whole family of LHV theories. In the past, 
many LHV models for Bell's tests with PDC have appeared in the literature 



13, uM. However, those models were mathematical constructs without too 



much physical content and, therefore, gave no hint as how to improve the 
experiments in order to make them reliable. In contrast, the LHV model 
which we present in this article gives a specific testable prediction: if that 
model (or a similar one) is true, necessarily there is a minimal light signal 
intensity which may be reliably detected. That minimal intensity depends on 
the geometry and characteristics of the optical devices (lasers, nonlinear crys- 
tals, lenses, etc.) and also on the properties of the photon counters used (in 
particular their quantum efficiency). The model disagrees with quantum me- 
chanics only for specific combinations of detection efficiency and experimental 
set up, and therefore determines the domain where a discrimination between 
quantum mechanics and local realism may be possible. Consequently, the 
model provides useful information about the domain where a reliable Bell's 
test, using PDC photon pairs, may be performed. Such a test will not rely 
on the plausibility of some additional assumptions, as has been the case in 
all experiments performed up till now. 
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In addition, our model refutes the wisdom that the quantum zeropoint 
field (the quantum vacuum fluctuations of the electromagnetic field) cannot 
be real because in that case it would saturate all detectors. Indeed, there is a 
long-standing controversy about whether the zeropoint field is real, or merely 
an artifact of the quantization procedure; during the 1940s these contrasting 
opinions were expressed by Casimir and Pauli respectively. If the fluctuations 
are not real it is difficult to understand phenomena like the Casimir effect or 
the Lamb shift. However, if they are considered real two big problems arise. 
One of these is the huge gravitational effect that it would produce at cosmic 
scales (more than 10 100 times greater than the known upper limits on the 
cosmological constant). We shall not be concerned with this problem here. 
The other difficulty is to understand how a so big radiation (about 10 5 w/cm 2 
in the visible spectrum alone) does not blind photon detectors, so making 
it impossible to detect single photon signals (as weak as 1 eV). Our model 
shows explicitly that the latter difficulty may be overcome, if we assume that 
actual photon counters have a detection threshold such that they detect only 
the radiation which surpasses it. The threshold may be considered as arising 
from the fact that atoms (and other quantum systems), in their ground state, 
are immersed in, and in equilibrium with, the zeropoint radiation, so that 
only radiation above the zeropoint level excites them. From this point of 
view, the quantum rule of using normal ordering in the calculation of photon 
absorption probabilities is seen as a mathematical procedure which takes 
account of the detection threshold. 



2 The Wigner function and the local hidden 
variables 

According to Bell, the crucial difference between quantum mechanics and 
LHV theories occurs in experiments where correlations between two particles 
at space-like separation are measured (Einstein, Podolsky and Rosen, or 
EPR, experiments [[HJ). Any LHV model should contain hidden variables 



A, with a probability distribution p(A), giving the following single and joint 
detection probabilities: 

Pl = J P (x)p 1 (x,4> 1 )dx, (i) 
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p 2 = I p(X)P 2 (\,<t> 2 )d\, (2) 

pia = J p(X)Px(X, 0i)P 2 (A, 2 )rfA. (3) 

where (f>i and 02 represent controllable parameters of the experimental setup 
and Pi (A, (pi), P 2 (A, <p 2 ) are some functions. If no further restrictions are put, 
a model resting upon Eqs.(|l|) to (|3p would be always possible. However, in 
order to have a LHV model, the functions Pi (A, 0i) and P 2 (A, (j> 2 ) should be 
probabilities and p(A) a probability distribution; consequently the following 
conditions should be fulfilled 

p(A) > ; J p(X)d\ = 1. (4) 

O<Pi(A,0i), P 2 (A,0 2 )<1. (5) 

In some particular instances a method to construct an explicit LHV model 
is to use the Wigner function. We consider a simple example. Let us assume 
that we have two particles initially (time t=0) in a state represented by the 
Wigner function W(xi,pi,x 2 ,p 2 ). Now the particles evolve freely until time t 
so that, according to the quantum rules, the Wigner function becomes W(x r 
Pit/m, Pi,x 2 — p 2 t/m, p 2 ) (remember that the quantum (Moyal) equation 
for the evolution of the Wigner function becomes the classical (Liouville) 
equation in the case of free particles). Finally we assume that at time t we 
perform local measurements on the particles such that we get the answer 
"yes" for the first (second) particle if it lies inside a region Ri(0i) (R 2 (4> 2 ))- 
We assume that the region Rj, (i — 1, 2 in the rest of the paper) is defined de- 
pending on one (or several) controllable parameters, (pi, of the experiment. In 
these conditions the probability of having the answer "yes" for both particles 
is 

P12 = I ,Pl,X2-— ,p 2 )©[xi G Pi(0i)]0[xi G P 2 (0 2 )]rf 3 xirf 3 x 2 , 

J m m 

(6) 

where 0[xj G Ri\ is the characteristic function of the region Pj, that is 
=1 (0 =0) if Xj belongs (does not belong) to the region P«. We see that the 
Wigner formalism provides an explicit LHV model for the experiment if the 
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Wigner function is positive. In fact, in this case equation (|6|) is a particular 
case of equation (0), the initial positions and momenta {xi, pi, x 2 , P2} play- 
ing the role of the hidden variables A. Explicit use of the Wigner function 
in order to obtain a LHV model was made by Bell himself [[UJ. In that 



paper, Bell showed that the situation proposed by Einstein, Podolsky and 



Rosen |15j had a positive Wigner function, and hence he was able to obtain 
a LHV model in the form outlined above. The problem is that only rarely 
the Wigner function of a two-particle state is positive, so this method cannot 
be generalized. 

If we pass from material particles to photons, then the Wigner function is 
quite frequently positive definite. In particular this is the case for the photon 
pairs produced in the PDC process ([Tj]] to JS2J). As it is well known, in 
the Wigner formalism the operators of creation, a k , and annihilation, a^, of 
photons become random complex amplitudes, a£ and a^, respectively. As the 
Wigner function is positive in this case, the amplitudes may be interpreted as 
that of a real random electromagnetic field, that is, these amplitudes become 
the local hidden variables, A, of our model with the Wigner function playing 
the role of the function p(A) entering in Eqs.(|l|) to In the case of photons, 
however, the detection probability has the form of a characteristic function, 
similar to 0(x« e R4) above, only if we use the Glauber function commonly 
labelled P({«k, but not if we use the Wigner function. Indeed, in 

quantum optics, the P function is assumed to correspond to the classical 
probability distribution of the radiation amplitudes. Therefore a LHV model 
parallel to that of Bell should use the function P, rather than W, and it is 
indeed well known that, whenever the P function is positive, the experiment 
does admit a classical (LHV) interpretation. 

If we insist on using the Wigner function, then it is not enough that it 
is nonnegative in order to get a LHV model. In addition we must introduce 
some non-negative functions of the random variables {a^a^K which play 
the role of the functions Pj(A, fa) of Eqs.([I]) to (Q) . Those functions should 
fulfill two conditions: a) They should be local, that is each function must de- 
pend only on the value of the electromagnetic radiation entering the detector 
during a detection time-window, b) They should give results so close to the 
quantum-optical predictions that the model is compatible with all performed 
experiments. Fulfilling these two conditions at the same time is not a trivial 
matter, and it is the goal of the present paper to present such a function. 
So most of the article will be devoted to exhibit a model of detection. Our 
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LHV model for the detection will rest upon the idea that the photodetector 
performs some average of the electromagnetic field during the detection win- 
dow and within the detector volume, and that a photocount will be produced 
when this averaged quantity exceeds a certain threshold value. 

In summary, our previous work ([[17]] to f22|) has shown that the Wigner 
function of the electromagnetic field at the detectors, derived from radiation 
produced by PDC, is positive for all performed experiments, and hence it 
provides an explicit LHV model for the production and propagation of this 
kind of radiation. The purpose of the present paper is to develop a LHV 
model for the detection, therefore completing an explicit LHV model for all 
PDC experiments. 

In section 3 we briefly describe the generation and propagation of PDC 
light in terms of the Wigner formalism. In section 4 we analyze the detection, 
clarifying why even with a positive Wigner function it is not trivial to get 
a LHV model. Section 5 is devoted to our detection model, the core of this 
paper. In section 6 the explicit calculation of the probability distribution of 
the LHVs of the model is done. In sections 7 and 8 we prove the agreement 
of this LHV model with the quantum mechanical predictions in the low 
efficiency limit. Finally in section 9 we discuss possible empirical tests of 
our model versus quantum optics. 



3 Wigner formulation of PDC: generation and 
propagation of light 

In parametric down-conversion (PDC) an UV laser pumps a non-linear (% 2 ) 
crystal and visible light is produced at a small angle (less than about 10°) 
respect to the incoming laser beam. The relevant fact is that there is a strong 
correlation between light beams produced at conjugate directions. In partic- 
ular if two detectors are placed in appropriate positions in order to detect the 
PDC light, coincidence counts are seen with a very high degree of temporal 
correlation (better than 10~ 13 s). The standard interpretation, resting upon 
the current quantum (Hilbert space) analysis of the phenomenon, is that 
each photon of the laser (u> , k ) splits into two conjugated ones, (u>i, ki) 
and (uj 2 , k 2 ), always satisfying the matching conditions: u> = u>i + uo 2 and 
ko ~ k\+ k 2 . Conjugate beams have some interesting coherence properties 
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which may be interpreted considering that the two photons of a pair are 
entangled. This is why PDC is considered a typically quantum phenomenon. 

As is well known the Wigner function formalism may be used for the 
study of quantum optics and it is completely equivalent to the more com- 
mon, Hilbert space formalism. Nevertheless the picture offered by the Wigner 
function is quite different from the standard one. It emphasizes the wave as- 
pects of light whilst the Hilbert space formalism stresses the particle aspects 
(photons are "created" in the source and "annihilated" at the detectors). 
A well known property of the Wigner function is that it exhibits explicitly 
the vacuum fluctuations of the electromagnetic field, also named zeropoint 
field (ZPF). In the Wigner formalism, once we take the ZPF as real, the 
phenomenon of PDC appears as a result of non-linear coupling, inside the 
crystal, between the laser beam and the ZPF. A picture emerges where every- 
thing looks like classical (Maxwell) wave optics, the Wigner function being 
the probability distribution of the field amplitudes. Also the electromagnetic 
radiation (including the ZPF) propagates according to Maxwell theory. This 
provides a straightforward wave picture of the generation and propagation of 
light for all PDC experiments (|[T/J to |22|]). Crucial for this interpretation is 



that the Wigner function is positive definite. In contrast the detection can- 
not be interpreted classically in the Wigner formalism, at least not trivially. 
The problem of detection will be considered in subsequent sections; here we 
shall summarize the main features of generation and propagation within the 
Wigner formalism. 

In the Heisenberg picture of the Hilbert space formalism, the way to 
compute single and joint detection probabilities (or rates) is based on the 
following expressions (modulo some constants related to the efficiency of the 
detectors): 

ptKiOlE-fatdE+fatiM, (7) 



p\ 2 cc (OlE-^tjE-ir^t^E+ir^t^E+^ttW). (8) 

(We shall use the superscript q when we refer to quantum mechanics and m 
for our LHV model). E(r{,ti) is the electric field operator at the position 
and time (r i; ti), and, as usual, we write E as a sum of positive and negative 
frequency parts, E + and E~ . These operators carry all the dynamical infor- 
mation of the system because in the Heisenberg picture the state remain fixed 
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whilst the operators evolve in time. This is the reason why these expectation 
values are computed with the initial (vacuum) state |0). 

In PDC, when we apply the Wigner function formalism, the creation and 
annihilation operators a k and contained in E~ and E + transform into 
random complex (c-number) variables a k and a k (k labels the normal modes 
of the electromagnetic field). Therefore, the operators E~ and E + transform 
into random functions E~ and E + (which we represent without hat) of these 
variables. What is the probability distribution for a k and «k? Because we 
are working in the Heisenberg picture it may be seen (Jl7| to f22|) that it is 



the Wigner function for the vacuum state |0), that is, the Gaussian function 

^(K,« k })=n^- 2|Qk|2 . (9) 

Then the electromagnetic field, now represented by the random functions E~ 
and E + propagates in a totally classical way when passing through all optical 
devices placed between the source and the detectors. 

The image of PDC that we get when studying the process with the Wigner 
formalism is similar to the classical one. Besides the laser, the crystal is 
pumped by the vacuum field, its positive frequency part now given by (for 
simplicity, we shall represent from now on the set of amplitudes {a k , a k } by 
a): 



ffr.i^l^E ^ a k e ikr -^. (10) 




This means that there is not only one input field on the crystal and therefore 
the PDC field is produced in a "classical" way as a consequence of the non- 
linear coupling between the laser and the ZPF. The PDC field has been 



computed in Refs. |21| and |22j , where it was shown the influence of the size 
of the crystal and the radius of the pumping on the matching conditions as 
well as the appearance of the typical rainbow in form of cones. A simplified 
expression of this field is 
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where three terms can be recognised. The first one is just the ZPF that 
crosses the crystal without any modification. The second one is the new 
field produced by the crystal as a consequence of the non-linear coupling 
between the laser and the ZPF. g is the coupling parameter which in general 
depends on u — u; for the sake of clarity we shall not write explicitly this 
dependence. The third term is similar to the first one; it consists again of 
ZPF, but a little modified. It is necessary to consider this term because 
the detection probability goes as g 2 . In order to make clear that we are 
dealing with a stochastic field, we have written in E {+) (r, t; a) the explicit 
dependence with the set of random amplitudes, a, which means that the PDC 
field has different values for different realizations of the ZPF. The amplitudes 
a obviously represent the hidden variables of our model. 
It is convenient to expand the PDC field in plane waves 

EM(r,t;a) = £ (^g^ ' & e~^\ (12) 
where the amplitudes /3k are linear functions of a, fulfilling the conditions 



IS : 



(APu) = Wif) = ; = o Mk'. (13) 

We end this section pointing out that, in the Wigner function formal- 
ism, entanglement appears as a correlation between the ZPF associated to 
two different light beams, whilst classical correlation just involves the signal 



leaving the ZPF untouched pOfl. This fact is becoming increasingly clear 



even out of the context of the Wigner formalism [23 



4 The Wigner formulation of PDC: detection 

In order to complete the Wigner function approach to PDC experiments, it 
is necessary to get the expression for the detection probability in terms of the 
electromagnetic radiation arriving at the detector. For the sake of clarity let 
us begin with the ideal case where the radiation field is represented by a single 
mode. In this case, working in the Heisenberg picture, the single detection 
probability is, according to quantum optics (for simplicity we shall omit the 
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subindex i and refer to the detector 1; obviously, the same expressions holds 
for the detector 2): 



Pi « ($\bl(r,t)b k (r,t)\$) = [bl(r,t)b k (r,t) + 6 k (r, t)&t(r, t) - l] |$> 

(14) 

= / [|/3 k (r,t;«)| 2 - 1] d« ^ <|/3 k | 2 - (15) 

where |<&) is the initial state of the radiation, W(a) the corresponding Wigner 
function, and b k (r,t) ( & k (r, t) ) the time-dependent creation (annihilation) 
operators, /3 k and being the corresponding amplitudes in the Wigner for- 
malism. The first equality derives from the use of the commutation relations 
and the second is the passage from symmetrical ordered operators to the 
Wigner representation. In the latter expression we have exhibited the de- 
pendence of the amplitude /3 k on position and time, and on the initial ampli- 
tudes, a. The simbol () w means the average of the quantity inside weighted 
with the Wigner function. In what follows we shall omit the subindex w 
In the general case, from (|7|) it is straightforward to get for a point-like 



detector 17 



p?(r,*)oc Jw(a)[I(r,t;a)-I ] da=(I-I ), (16) 

where J(r, t; a) = c£qE + (t, t; a)E~(r, t; a) is the intensity for a realization 
of the field ( |12"D at the position and time (r, £), and W(a) is the Wigner 
function of the initial state. N is the number of modes (we should take the 
limit N — >• oo at some appropriate moment). Iq is the mean intensity of the 
ZPF, so that Eq . (|16D might be interpreted as stating that the detector has a 
threshold so that it only detects the part of the field that is above the average 
ZPF. Two remarks are in order: a) The intensity I contains a (possibly 
complicated) dependence on the initial amplitudes of all radiation modes, but 
Iq is a constant (compare with (|T5|)). We might understand this fact by saying 
that the detector "knows" the radiation actually arriving at a given time, that 
is signal plus noise (ZPF), but it is not a trivial matter to remove the noise. 
The quantum rule (|l^) is just to subtract the mean, a formal procedure 
which cannot be physical because it gives rise to "negative probabilities", 
as we shall discuss below, b) Strictly the integral in (|i~6D should involve 
all radiation modes; therefore some cut-off frequency is required in order to 
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avoid divergences. Nevertheless, most of the radiation modes are usually 
not "activated" (in usual quantum language we say that they contain no 
photons) and therefore they may be ignored. That is, the contribution of 
these modes to the average (I) equals the contribution to J so that ignoring 
them does not change the difference (I — Iq). We shall call "relevant modes" 
those which cannot be ignored. 

In the same way it is also possible to obtain the coincidence detection 
probability for two detectors, placed at (ri,£i) and (^,£2). We get from (||) 

pf 2 (ri, ti; r 2 , £ 2 ) oc / W(a) [I (r 1 , £1; a) - I w ] [I (r 2 , £ 2 ; a) - ho] da 

= ((/ 1 -/ 10 )(/ 2 -/ 2 o)), (17) 

where we see that the coincidence detection probability is obtained by mul- 
tiplying the intensities that arrive at the two detectors, after subtracting the 
threshold, and averaging this quantity with (§) for all the possible realiza- 
tions of the field. Writing I\q 7^ /20 we emphasize that the thresholds of 
both detectors may be different (in particular the relevant modes are usually 
different in the two cases). 

The relevant question for us is whether expressions ([16]) and (|17|) are 
suitable for a local realist interpretation. In fact, those expressions look like 
Eqs.(P to (^|) with W(a) playing the role of p(A) and / — I the role of P. 
However, I — I does not fulfil Eq.(^); in particular it is not always positive 
and therefore the quantum theory of detection prevents to get a trivial LHV 
theory from the Wigner representation. The problem is not the huge value 
of the zeropoint energy because the threshold intensity Iq cancels precisely 
that intensity. The problem lies in the fluctuation of the intensity. For the 
weak light signals of the experiments the fluctuations of / may be such that 
K/o. 

In order to study whether this important problem may be solved we shall 
begin showing that the removal of some idealities involved in Eqs. ([UJ) or 
( |TTD alleviates the situation. Firstly these equations were derived including 
only modes corresponding to a beam of almost parallel wave vectors. If this 
is not the case we should write Eq. ( |P7|) using the Poynting vector rather 
than the intensity. The direction of Poynting vector of the signal is well 
defined whilst that of the noise (the ZPF) is random with zero mean, which 
may make easier the discrimination. More important is the fact that the 
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detection probability should not depend on the instantaneous intensity at a 
point. In fact, it is natural to assume that the detection probability depends 
on all the incoming radiation entering the detector during the detection time 
window. We shall write Eqs. ([16]) and flT7|) in a more realistic form, with a 
different meaning for the average represented by (). We have 

pi = J W(a)Q 1 (a,<l>i)da , (18) 
P12 = J W{a)Qx{a, 0i)Q 2 (a, <f) 2 )da , (19) 

where 



Qi{a, fa) = -T- 1 - / dU / d 2 ri[Ii(a,<f)i,ri,ti) - I i0 ] = h(a, 
hvi Jo J a hvi L 



'H - HO 



(20) 

where Ii(a,(pi) — = Ii s (a,4>i) is the result of integrating the difference 
between the actual intensity and the intensity of the zeropoint field over the 
time window, T, and the surface aperture of the detector, A. We have divided 
by the typical energy of one "photon" so that Qi becomes dimensionless. In 
this way rji is the quantum efficiency of the detector. 

The relevant question is whether ([TJJ) and ([[£]) may now be considered 
a particular case of ([l]) to (§). If the answer is affirmative (negative) the 
formalism provides (does not provide) an explicit LHV model for the experi- 
ment. Actually the answer is not yet affirmative because we cannot guarantee 
the positivity of Qi. (For the additional requirement, Qi < 1, see below Eq. 
(p8[) , but this condition certainly holds for the low collection efficiencies of the 
experiments performed up till now). The problem is now alleviated by the 
time and space integrations in (|20|). Indeed, the fluctuations of the intensity 
are strongly reduced by averaging over space-time regions, as the Heisen- 
berg (uncertainty) relations show. But we can guarantee the positivity of Qi 
only in the limit of infinitely wide time-windows or infinitely large apertures, 
which is non-physical. Consequently we conclude that it is not possible to 
interpret directly the Wigner-function formalism as a LHV model for the 
PDC experiments. 

We shall devote the rest of the paper to describe a plausible model of 
detector that "works" in a strictly local way. 
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The detection model 



3 will now proceed to show the basic points of our model: 

1. We shall assume that the detector is formed by a set of individual 
photodetector elements, Dj, each characterized by a central frequency 
Uj, and a wave vector kj [ujj = |kj|/c), to which Dj responds. We shall 
consider the direction of kj to be normal to the surface of detector, 
which is taken as a cylinder of area A = nR 2 and length L. 

2. If the light beam contains frequencies into the interval (kVin.^max), then 
ujj G (^ miri) u; max ). Also, given two detector elements with frecuencies Uj 
and Ui, the following relation holds 



\u)j-ui\>^. (21) 

We shall suposse that the relevant quantity for the detection is not 
directly the electromagnetic field E^ + '(r, t;oc), but a filtered field (by 
Fourier transform over the detector volume and time window of de- 
tection). We shall define the filtered field corresponding to detector 
element Dj, as 

E { +\a) = J y dV £(+)(r,t;a)e-^ r +^. (22) 

This equation shows that will depend only on the radiation that 
crosses the detector during a time window, which is the locality con- 
straint on the model. 

By substituting Eq. flT2"|) into Eq. (|2"2"| ) and taking the origin of the 
reference system OXYZ at the center of the crystal with the OZ— axis 
as its axis, we obtain, after some easy algebra 



1/2 



sine 



T 



(23) 

The presence of the sine factor in Eq. ( p3|) implies that, for a given value 
of Uj, only frequencies whithin a range of width 2n/T centered at ujj 
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will contribute to the sum in (|23|), i.e. the photodetector element Dj is 
sensitive to radiation with frequencies in the interval (uj — ojj + 4p) 
with Au ~ 2tt/T. 

4. The maximun number of independent detecting elements is (see Eq. 
&) 

N^ — ^-, (24) 
Au T 

where Sou = c<j max — uj min pa 2tt/t, t being the coherence time of beam. 
Typical values for T (pa 1CT 8 s) and r (pa 1CT 12 s) give rise to N = 10 4 . 

Note that the modes which contribute to Ej are different to those cor- 
responding to , which implies that E^ and e\, ^ are uncorrected 
for jV k (see Eq.(|I3|)) : 

(E^s^) « o ; ( 25 ) 

This expression will be useful later on. 

5. Now we shall define the effective intensity, I(ct), which is obtained from 
the filtered fields ( |22"D in the form 

N 

7(a) = ceo^W+\a)Wf\a), (26) 

3=1 

and replace Eq.(EDI) by the expression 



P(a, 4>) = ((I(a, 0) - /o)6[/(a, 0) - I m ], (27) 

where Jo is the average of / for the ZPF. I m is some threshold intensity 
related to the "voltaje biass" of the detector, fulfilling the condition 
I m > Jo, and 9(x) is the Heavside function O(x) = 1 if x > 0, other- 
wise. We have considered all the complex dependence of the intensity 
on the initial amplitudes, a, and the controllable parameters in the 
experiment, 0. 

We see that taking ( — > r/TA/hu, our Eq.([T7|) becomes the standard 
quantum detection probability Eq.(Tl6|) if we consider instantaneous 
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detection (i.e. T — > 0) by a point-like detector (i.e. R — > 0, L — > 0) 
and remove the Heavside function. We would like to interpret Eq.fl27]) 
as the probability for a given realization of the field and it is analogous 
to Pi(A, 4>i) or P 2 (A, 2 ) in Eqs. (|l|) and (|j) for a value of the hidden 
variable. However, although the Heavside function ensures that the 
quantity P(a, 0) in (p7|) is positive, it does not fulfil the condition 
P(a, 0) < 1 in general. Therefore we should propose an expression 
which is positive, lower than 1, and reduces to Eq. (p7[) in the case of 
(\I(a,(fi) — I \ < 1. A simple expression achieving these goals is 

P(a,0) = (1 -e-« 7 ^- J °))e[7(a,0) - I m ], (28) 
which completes the definition of our model. 
6. The single detection probability predicted by the model is 



p™ = J W(a)P(a,(f>)da, (29) 
which can be expresed in the following equivalent form: 

p? = J p(l)P(l) dl , (30) 



where 

N 



p(!) = [ W(a)6[I -c€oY,E ( j h \a,0)E ( j \a,<f>)]da. (31) 
J 3=1 

Similarly we have, for the joint detection probability, 

P12 = J Pi 2 (7i,7 2 )Pi(7i)P 2 (7 2 )rf7irfJ 2 , (32) 

where 



Pn{h, h) = I W(a)5[h - ce Ef =1 E^\a, ^^(a, 00] 
x5[7 2 - ce Ef=iE%\a, h)E { ~\a, </> 2 )]da. 



(33) 



Eqs. (^2) and (^2|) are very convenient for the comparison between our 
model and quantum optics. As said above, the amplitudes a are the hidden 
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variables, W(a) playing the role of the function p(X) in Bell's formulation. 
The analogous of the probabilities P(X, 0) are 



P(a,<f>) = f P(/)5[/-ce o f;^ +) (a,0)4 \a,<j>)) dl. (34) 



The (complicated) dependence of the filtered field on the vacuum amplitudes 
should be derived from the quantum Wigner formalism for every experiment. 



This was made in Refs. |T7j to 

Now we shall study whether the predictions of our LHV model are com- 
patible with the results of performed experiments and to which extent they 
agree with the quantum predictions. For that purpose the following steps 
will be taken: i) To obtain the probability distribution, po(I), of the effective 
intensity of the ZPF in absence of any further electromagnetic radiation, ii) 
To obtain the probability distribution, p(I), for the radiation that arrives 
at the detector when it is illuminatated with the light coming from a PDC 
process. Hi) To compute single detection probabilities and iv) joint detection 
probabilities. 



6 Calculation of po(I) and p(I) 

The filtered field E + is Gaussian because it derives from a Gaussian ran- 
dom field, Ej~, under linear transformations. Consequently E + Ej ^ is an 
exponential random variable. Nevertheless, by virtue of / being the sum of a 
large number of independent random variables (~ 10 4 ), a version of the Cen- 
tral Limit Theorem applies, so that / is Gaussian to a good approximation. 
Consequently the full probability distribution is determined by the mean and 
the standard deviation. 

For instance, the probability distribution for the effective intensity when 
only the zeropoint field is present is 

_ 1 a-Tp) 2 

p (I) = -^e^T. (35) 
\Z2irao 

Let us calculate I and cr . By using Eqs. ( p3|) (making (3-^ — > a^) and 
(|26|), we have 
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7o = E J o„ (36) 

where 



hj 



h cAJl(k r R) 2 



L 



10 A 



sine 



(37) 



Here we have made use of the relation (a^a^,) = 5^/2, which can be easily 
derived from Eq. @. 

Taking the continuous limit Sk/-^o — > I d 3 k/(2ii) 3 , and changing to 



spherical polar coordinates (k X} k y , k z 



LO 00 00 

k z = — cos6> ; k x = — sen6 cosip ; k v = — sen.6 sem/>, 
c c c 



we get 



'- Oj 



2nch 
2(2vr) 3 



a; 2 da; 
c 2 c 



■sine 



x 



t/2 4J 1 2 (^sinfl J R) „ 

g?0 sinfl t~~ . c — r-r sine 



— (a;cos6 l — a;^) 



(38) 



'-smOR) 2 

The typical value of Too, uj being a frequency in the visible range, is 10 7 . For 
this reason we can approximate the first sine factor in the above expression 
by a delta function. After performing the integration in a;, we obtain: 



hj 



K 3 



Ac 2 Tn Jo 



tt/2 



ad sinfl 77^—7 —7- sine 



^(C08*-1) 



^sinfli?) 2 

Now, making the substitution v = Lujj(l — cos 9) /2c, we have 

4J 2 



(39) 



hj 



H 2 



2ttcTL Jo 



Lojj/(2c) 



dvsmc (v) 



1- 



(40) 



In this expression the upper limit of the integral is much larger than 1, 
because L ^> A, A being the typical wavelenght of the radiation. For this 
value of v the sine 2 factor is negligable. Thus we can safely extend the 
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upper limit to infinity. Also, the fraction AJ 2 (x)/x 2 can be approximated by 
1 m case that x 2 /A < 1. It can be shown that this implies the condition 



R < \j AL/87T 2 , a situation that stands for usual small detector radius. Under 
these conditions, we obtain 

hi = — ■ (41) 

03 AcTL v ; 

From (53), the total mean intensity is 



f^AcTL f^ AcTL 



Passing to the continuous, we finally obtain 



where u = (u max + u min )/2. 

Let us now proceed to the calculation of a . We have 



(42) 



T p« foul? _ huJ 2 5uj 
~ 2^ L min AcTL ~ 8ncL ' ^ ' 



-o 2 = a 2 )-7o = c 2 ^EE 



\^oj ^oj ^o* ^% / - (-^oj E 0j I \ h m h oi ) 



(44) 

Taking into account that we are dealing with Gaussian processes and that 
(E^E^) = (E^E^) = (V j, I), and (E^E^) « Oif ^ j (this follows 
from the fact that we are dealing with a set of iV practically independent 
stochastic variables), we find 



(45) 



Hence we obtain 



<T Q = <r 0j y/N=(I )yJ^, (46) 



and the final expression for po{I) is 
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where Iq is given in Eq. ([43D . 

Let us proceed to the calculation of the probability distribution for the 
effective intensity when there is a PDC signal present, p{I)- Such distribution, 
as that of po(I), is also a Gaussian. The argument is parallel to that leading 
to Eq. (f4TD starting from the fact that the filtered field (see Eq. (^)) is 
Gaussian and its statistical properties are close to those of the ZPF, except 
for the greater intensity of the former for some modes of the radiation. The 
Gaussian character of the PDC radiation has been discussed in detail in our 
previous papers (jT7| to [22]). But there is also another argument derived from 



the fact that the counting statistics corresponds to a Gaussian (Glauber) 
P function (see Jj4j] ). Now, it is known that the Wigner function is the 
convolution of the P function with the Wigner function of the vacuum, which 
leads to a final Gaussian Wigner function. The mentioned convolution has 
an interesting physical interpretation: It corresponds to a random variable 
which is the sum of two uncorrelated random variables. It is as if the signal 
is superimposed incoherently to the ZPF when P> 0, a situation usually 
named "classical light". 

The relation between the mean and the standard deviation is completely 
analogous to that of the ZPF alone, and the distribution may be written 

_ \— (T-<T» 2 t \/~ _iZzIozik£l 

p(I) = JU_e _ ^ = -JLW" -oV (48) 

27Ti y/ Z7T Iq 

where we have defined the signal mean effective intensity by 

7.= <J)-7 , (49) 
and taken into account that I s -C Iq in the denominator of the exponent in 



([48D . This can be easily demonstrated for typical values of L, r, and A in Eq. 



(EST), and by considering the intensity arriving at the detection system in a 



usual PDC experiment. 
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7 The single detection probability 



We want to start by comparing the prediction of our model with that of 
quantum optics. The quantum optical prediction for the single detection 



probability is given by Eq. ([181). By using the properties of the delta function 



this expression can be written in a different way: 

p\ = JL r dl (I - J ) / da W(a)5[I(a, 0) - I] 
hu Jo J 



where 



p(I) = J da W(a)8[I(a, <f>) - 1\. (51) 

Now let us analyze the predictions of our model. They are given by 
Eq. (^) where P(I) was defined in Eq. (|28|) and p(I) in Eq.([H]), the latter 
calculated in section 6. After some easy algebra, we obtain: 



In the above expression the important parameters I s /(Tq, (I m — Iq)/oqi 
and I s ( enter. Let us now consider the two following situations: i) (I s <C 1 
and (<Jq <C 1 (linear approximation), and ii) (I s ^> 1. 

i) (I s <C 1 and (ao <C 1. This should be the normal situation in ex- 
perimental practice. In order to compute the detection probability in this 
case, we take into account that the relevant values of J — Jo in the inte- 
gral of Eq. (|30"D are close to I s , which allows to make the approximation 
1 - e-« 7 - 7 °) C (7 - 7o) ■ We obtain 

p? = ^erfc[ /m ~ /o " /s ] + -^e" (Im "-r )2 . (53) 
Now, we may choose I m so that 

7 S - (J m - 7 ) > (T , (54) 
With these two approximations we arrive to the following result: 
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P? = d s . (55) 

The important point is that this result is very similar to the quantum one, 
the only difference being the presence of the effective intensity of the signal, 
I s , instead of the actual intensity, (I S )/AT, in the quantum formula. In 
experimental practice a lens is placed in front of the detector in such a way 
that the signal field has spatial coherence on the surface of the lens. The 
condition for having spatial coherence is 

d\ > RiR c , (56) 

d being the typical distance between the nonlinear medium (with an active 
radius Rc) and the detector; Ri is the radius of the lens. It can be demon- 
strated (see Appendix) that this property of the signal field implies that 

/. = f. (57) 
Hence, by substituting Eq. fl5"T| ) into Eq. fl55|), we have 

P? = d s = ^-(l), 
hv 



a result that coincides with Eq. (|18|) . 

If is = strictly, in sharp contrast with quantum optics, our model 
predicts the existence of some counts in any detector even in the absence of 
signal. In fact, from Eq. ( |53|) we get 



(a 



2 



ri£r* = -7==e 2CT o , (58) 

V 27T 



which is very small if the condition 

i m -i >a , (59) 

is fulfilled. In this case there is no conflict with experiments, because we may 
interpret the counts without signal as a part of the dark rate of the detector. 

ii) (I s ^> 1. By taking this limit in Eq. (^) one can easily check that 
under this condition p™ = 1 (provided that the constraint ( |54|) holds true). 
That is, the detector saturates when the intensity is very high, and gives a 
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count in every time window, a fact not in disagreement with experiments (or 
with quantum-optical predictions). 

However, the remarkable agreement with quantum-optical predictions 
is not enough to guarantee that we have arrived to a LHV model for the 
detection probability. This is because in our model there are two constraints 
that must be fulfilled: the constraint (0), required for the linearity of the 
response at low efficiency, Eq.(|55]), and (|59D , needed for the smallness of the 
dark counting probability, Eq. (|58|) . For these two conditions to be fulfilled it 
is necessary that I s ^> <7o, a condition which makes difficult to construct a 
LHV model for detection and, at the same time, what gives predictive power 
to the model, as will be discussed in detail in the final section. 



8 The joint detection probability 

The quantum mechanical prediction for the joint detection probability is 
given by Eq. (|1^). As before, using the properties of the delta function this 
expression can be written in the following way: 



P Q i2 = Jo°° dh Jo°° d~h I da W(a)5[h(a, fa) - h]5[I 2 (a, fa) - I 2 ] 
xft(Wio)£(W*) 

= r °di x r , d/ 2 p(/ij 2 )^(Ji-/ w )^(/ a -/ 20 )= 1 ^</ 1 ./ 2 .>, (60) 

jo Jo hv\ hv 2 n z viv 2 

where 

p(h, I 2 ) = J da W(a)8[h(a, fa) - h]5[I 2 (a, fa) - I 2 \. (61) 

Now let us analyze the predictions of our model for the joint detection 
probability. They are given by Eq. ( |32|) where Pi(Ii) was defined in Eq. (f2lf) 
and p(hj 2 ) in Eq. ©: 

poo roc m _ — — — — „ ™ 

P?2 = / / Pu(h, h){l ~ e-^-^l - e-^-^dhdh. 

J Ilm 1 2m 

Because p(I\, I 2 ) is a doble Gaussian function, it is defined by the mean 
values of its marginals, their standard deviations and the correlation function 
((I i - I u - Ii )(I 2 - l 2s - 7 2 o))- In the case that o x = a 2 = cr , then 
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p(h,I 2 ) = (2ira 2 Q 



2\-l 



-1/2 



exp 



/2 + /|_2/ 1 / 2 (/ 1 / 2 )/ ( T 2 



(62) 

where we have defined l\ = l\ — l\ s — Iio and J 2 = 1 2 — 1 2s — ^20 • 

Now, by considering the limits QI is <C 1 and I is — (I im — I i0 ) ^> cr , we 
arrive to the following result: 



pS ~ / °° d/i / °° di 2 p{h,i 2 )Uh -hoMh - ho) 

= ClC 2 ((/l-/20)(/l-/20)), 

(63) 

a result very similar to the quantum prediction, Eq. fll~9|) . 

The parameters in p(I\, 1 2) are easily obtained from the marginals and the 
correlation function by making equal their values to the quantum mechanical 
predictions: 



00 roc 



dl 2 / I lP (h,I 2 )dI 1 
JO 

/■ :x. 

dl 2 I 2P (I U I 2 )dh 



dhp{h)h = I w + I u , (64) 
dl 2 p{l l )l 2 =l w + l 2s) (65) 



poo poo „ _ _ „ _ _ 

/ dl 2 dhhl 2 p(l u h) = (hl 2 ). 

Jo Jo 



(66) 



9 Discussion: Empirical tests of the model 

The predictions of our LHV model agree with those of quantum theory for 
all PDC experiments with low efficiency, as stated above. Therefore our 
model violates all "Bell's inequalities" empirically tested up to now, those 
inequalities having been derived from local realism plus auxiliary assump- 
tions. Consequently the model violates the auxiliary assumptions, although 
we will not discuss this point here in detail. At higher efficiency, the model 
departs from the conventional quantum theory in that it predicts a nonlinear 
response of the detectors (see eqs.(|30|) and (^)), or a high dark rate (see 
ec L-(il))) or both. This is the feature that prevents the violation of a genuine 
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Bell's inequality (that is, one involving no auxiliary assumptions in addition 
to local realism). 

An obvious disproof of our LHV model would be achieved if a PDC ex- 
periment violated a (genuine) Bell's inequality; this is what is usually called 
a "loophole-free" test. But, at its most optimistic, such a test, using para- 
metric down-converted photon pairs, lies in the remote future. In any event, 
such a test would necessarily involve measurements of both single and co- 
incidence counts |1(J, and should avoid any background subtraction. The 



latter condition derives from the fact that our model predicts the existence 
of a fundamental dark rate in photon detectors and there is no reason why 
it should satisfy a Bell's inequality if that rate is subtracted. 

But simpler tests may be found for our model, since it has a predictive 
power greater than that of quantum optics concerning the behaviour of de- 
tectors. Nobody would claim that quantum mechanics had been violated if it 
were discovered that there is a dark rate, or that the response of the detector 
to the signal is non-linear, even though both of these are in disagreement 
with the quantum theory of measurement. These facts would be attributed 
to imperfect functioning of the detectors, and the imperfections would be 
considered as technical problems of no relevance to the testing of quantum 
theory. In sharp contrast, our model establishes rather stringent fundamental 
constraints on the functioning of detectors, precisely because the zeropoint 
field is taken as real. And the reality of the ZPF is an unavoidable con- 
sequence of taking the Wigner function as the probability distribution of 
hidden variables, which is the central idea of our model. 

The constraints posed by our model were already mentioned in section 7, 
and they arise from the necessity of making compatible the conditions ( pUD 
and ([54]), which imply 

7 S » (To- (67) 
Hence, if we take into account (pESf) and ( pES| ) we get 

- hcj 2 , , 

» W5T (68) 



This means that, in our detection model, there is a minimal effective intensity 
of the signal which may be reliably detected, a constraint absent in the 
quantum theory of detection. This is the constraint which may be put to 



25 



empirical test (remember that this constraint is required only if we demand 
a low dark rate). Now we shall analyze the consequences of the constraint. 

As a consequence of the fact that the signal field has spatial coherence 
on the surface of the lens, the intensity of the incident signal is amplified by 
a factor 

b 2 = 7r 2 Rf /A 2 / 2 , 

/ being the focal distance. On the other hand, the zeropoint field is not 
modified by the lens, which is evident because of the fact that energy cannot 
be extracted from the vacuum. 

84 % (91 %) of the total intensity is concentrated within the first (second) 
ring of the difraction pattern, with a radius R = a x (f\/2Ri) = a\/A r , 
A r = 2R[/ f being the relative aperture of the lens, and a = 1.22 (2.23) for 
the first (second ring) f25[. Consequently, the optimus radius of the detector 
is given by R. After that we may write constraint (^) in terms of the 
intensity, I/at, arriving at the aperture of the detection system as follows 



- 2T m 2 \ 2 f 2 

(in = b~ 2 I s » 



Air 2 RfcL\rrT 



Still, we may write the constraint in terms of the most direct empirical quan- 
tity, namely the single counting rate, and we get 

Rate = ^f IN » ? f * (69) 

where rj is the quantum efficiency of the detector and we have used the 
equality uJ = 2ttc/X. It may appear that the detection rate could be as low 
as we want by just increasing the radius, R;, of the lens, but this is not true 
because the condition ([56]) should be also fulfilled. If we combine this with 
)BT) we obtain 

Rate » ^f 2R c (7Q) 



The constraints ( p9|) and ( |70D , which put a lower bound to the single rate 
which may be used in reliable experiments, is the most dramatic prediction of 
our model. The existence of a lower bound to the single rate for the reliability 
of PDC experiments cannot be derived from (conventional) quantum theory. 
If we put typical parameters of the detector, that is rj ~ 0.1, Ri, L and 
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f of order of fractions of a centimeter and a time window T pa 10 ns, and 
use typical values of the wavelength, A ~ 700 nm, and bandwidth, AA pa 10 
nm (which gives a coherence time r pa 1 4 ps) we get a minimal counting 
rate of the order of 10 5 — 10 6 counts per second. This figure is not far from 
the one appearing in actual experiments. In any case, the model requires 
improvements in order to be able to make more accurate predictions, a work 
which is in progress. 



10 Appendix 



In order to demonstrate Eq. fl57p we shall first substitute Eq. fl2"2]) into Eq. 
(|26|). We shall express the electric field as the sum of two terms (see Eq. 



EM(r,t;a) = E^ + \r,t;a) + Ei + \r,t; 



a 



(71) 



where E^ + \r, t; a) is the part of the field which is superimposed to the zero- 
point field Eq (r, t; a). The calculation of the effective intensity correspond- 



ing to the filtered zeropoint field has been computed in Eq. (f43|). Let us now 
focus on the calculation of I s , which is given by the following expression: 



= G^fer) 2 '^olv^dV Jy^dV'ff dt / T dt>(E(+\v,t;a)E(-\r',t>;a)) 
x £\ Q- i ~i-( z - z ')+ icj j( t - t ') 

(72) 

The spatial coherence implies that the integrals over the surface of the de- 
tector are equal to irR 2 . Passing to the continuous in the summation in j 
Eq. ( [72] ) transforms into 



I s = (±) 2 csq Cl'il dz £$dz? So 'dt £ dt'(E^(z,t;a)Ei-\z',t';a)) 



27T ■'k'mm 



-i-t(z-z')+iwj(t-t') 



^-Bhh S-l/2 dz dz' Jo dt / T dt'(EM t ; a)E^ (z>, t'; a)) 



-)]. 



(73) 
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Now, we make the substitution sinc[^(t -If- '^f)\ -> ^5(t - 1' - ^f), 
and perform one integration on time. We have 

7, = ^I^dzJ^dz'J^dt(E^Kz,t-a)Ei-\z',t-^)) 

= ^f^dztfdt(Ei+)(z,t;a)E^(z,t;a)). (74) 
& = t\ I-l/2 dztf dt ce (Ei + \z, t; a)E^(z, t; a)). 



This result coincides with the quantum one, Eq. (18), when the spatial 
coherence is taken into account in that equation. 
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